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We successfully carried out the first high-altitude balloon flight of a wide-field hard X-ray coded-aperture telescope 
ProtoEXISTl, which was launched from the Columbia Scientific Balloon Facility at Ft. Sumner, New Mexico on 
October 9, 2009. ProtoEXISTl is the first implementation of an advanced CdZnTe (CZT) imaging detector in our 
ongoing program to establish the technology required for next generation wide-field hard X-ray telescopes such as 

f^v the High Energy Telescope (HET) in the Energetic X-ray Imaging Survey Telescope (EXIST). The CZT detector 

plane in ProtoEXISTl consists of an 8 x 8 array of closely tiled 2 cm x 2 cm x 0.5 cm thick pixellated CZT crystals, 
each with 8x8 pixels, mounted on a set of readout electronics boards and covering a 256 cm 2 active area with 

i^^ 2.5 mm pixels. A tungsten mask, mounted at 90 cm above the detector provides shadowgrams of X-ray sources in 

the 30 - 600 keV band for imaging, allowing a fully coded field of view of 9° x 9° (and 19° x 19° for 50% coding 
fraction) with an angular resolution of 20'. In order to reduce the background radiation, the detector is surrounded 
by semi-graded (Pb/Sn/Cu) passive shields on the four sides all the way to the mask. On the back side, a 26 cm x 
26 cm x 2 cm CsI(Na) active shield provides signals to tag charged particle induced events as well as > 100 keV 

. ^h background photons from below. The flight duration was only about 7.5 hours due to strong winds (60 knots) at 

float altitude (38-39 km). Throughout the flight, the CZT detector performed excellently. The telescope observed 
Cyg X-l, a bright black hole binary system, for ~ 1 hour at the end of the flight. Despite a few problems with the 
pointing and aspect systems that caused the telescope to track about 6.4 deg off the target, the analysis of the Cyg 
X-l data revealed an X-ray source at 7.2a in the 30-100 keV energy band at the expected location from the optical 
images taken by the onboard daytime star camera. The success of this first flight is very encouraging for the future 
development of the advanced CZT imaging detectors (ProtoEXIST2, with 0.6 mm pixels), which will take advantage 
of the modularization architecture employed in ProtoEXISTl. 
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1. Introduction 

The high energy sky is active and dynamic, full 
of inherently variable sources as well as extremely 
luminous Gamma-ray Bursts (GRBs). Wide-field 
hard X-ray imaging is perhaps the only effective 
way to capture these unpredictable, energetic ce- 
lestial events. The current GRB mission, Swift, to 
study GRBS, employs a large area (0.5 m 2 ) array 
of CdZnTe (CZT) detectors, each 4 mm x 4 mm x 
2 mm as separate pixels for coded-aperture imag- 
ing. The previously proposed EXIST telescope 
would have had 4.5 m 2 with 0.6 mm pixels [2]. A 
smaller (1024 cm 2 ) CZT imager incorporating the 
same 0.6mm pixels has now been proposed for the 
MIRAX experiment on the Lattes mission to be 
launched by Brazil [3] . The combination of a large 
area and fine pixels imposes serious challenges in de- 
tector design, packaging, and operations, given the 
limited power and space available in space missions. 

We have been pushing the technology of advanced 
CZT imaging detectors required for future genera- 
tion of wide-field hard X-ray imagers through our 
ongoing, three-part, balloon-borne wide-field hard 
X-ray telescope experiment, ProtoEXIST. We have 
successfully carried out a first high-altitude balloon 
flight of ProtoEXISTl, the first implementation of 
the program, which was launched from the NASA 
Columbia Scientific Balloon Facility (CSBF) at Ft. 
Sumner, New Mexico on October 9, 2009. ProtoEX- 
IST1 is a wide-field hard X-ray coded-aperture tele- 
scope, built to develop an efficient modularization 
and packaging architecture for a large area fine pixel 
CZT imaging detector (256 cm 2 with 2.5 mm pix- 
els) and to demonstrate its performance in a near 
space environment. The X-ray telescope including 
the CZT detector performed very well during the 
flight. Here we report the detailed performance of 
the X-ray telescope and the imaging CZT detector 
during the flight. Full details of the detector and tele- 
scope system, including the assembly of the coded 
mask, the shields, the onboard calibration source 
and the interface to gondola will be presented in sep- 
arate papers [4,5]. 

2. ProtoEXISTl Instrument 

ProtoEXISTl consists of a CZT detector plane, a 
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Fig. 1. Schematic view of ProtoEXISTl. The mask and the 
upper side shields are mounted on the top cover of the pres- 
sure vessel (PV), and the rest of the system is mounted in- 
side the PV. See also Table 1. 

Tungsten mask, passive side shields and a Csl rear 
active shield. The schematic view of the ProtoEX- 
ISTl telescope is shown in Fig 1. The key parame- 
ters of ProtoEXISTl are summarized in Table 1. 



2.1. CZT detector plane 

The CZT detector in ProtoEXISTl is a first im- 
plementation of a close-tiled fine pixel imaging CZT 
array in our program. Fig. 2 shows the detector as- 
sembly and the fully integrated detector plane in 
the flight configuration. An initial laboratory perfor- 
mance and the assembly procedure of the CZT de- 
tector plane for ProtoEXISTl are reported in [6] in 
detail. Here we review the modularization and pack- 
aging architecture and operational modes briefly. 

The full detector plane covers 16 x 16 cm 2 with 8 
x 8 CZT crystals, providing 64 x 64 imaging pixels 
with 2.5 mm pitch. All of the 64 CZT crystals are 
tiled with a uniform gap (~ 900 /im for 19.5 mm x 
19.5 mm CZT), allowing the total active area of 15.6 
x 15.6 cm 2 . This is a factor of > 8 times larger in 
area than any previous tiled-pixellated CZT imagers 
of which we are aware. 

The basic building block is a Detector Crystal 
Unit (DCU) that consists of a 1.95 cm x 1.95 cm x 
0.5 cm CZT crystal with 8x8 pixels bonded to an 
interposer board with a RadNET ASIC [7] under- 
neathfor readout (Fig. 2a). The traces on an inter- 
poser board connect the 2-D array of 8 x 8 anode 



pixel pads on a CZT to a 1-D array of 64 channel 
input pads of the RadNET ASIC. In the next stage 
of assembly, a 2 x 4 array of DCUs are mounted on 
two vertically stacked electronics boards to make 
a Detector Crystal Array (DC A) (Fig. 2b). The 
DCA electronics boards contain four Analog-to- 
Digital Convertors (ADCs, AD7685BRM) and a 
Field Programmable Gate Array (FPGA, Altera 
Cyclone II EP2C20F256C7). The FPGA is pro- 
grammed to have four independent data processing 
channels, each of which handles two DCUs with 
one shared ADC. We closely tile 2x4 DC As on a 
motherboard, the FPGA-Controller Board (FCB), 
to complete the detector plane or Detector Module 
(DM; Fig. 2c). The FCB contains another larger 
FPGA (Altera Cyclone II EP2C20F484C7) to pro- 
cess 32 independent data streams (1 per 2 DCUs) 
from the FPGAs of the 8 DCAs. A network card 
(NetBurner 1 , MOD5282) is mounted on the FCB 
to collect the science data from the FPGA and 
other HouseKeeping (HK) information, and pass 
them to a flight computer through an Ethernet port 
(see §2.3 and [5]). The FCB also reserves 8 TTL in- 
put lines for event tagging with 241 Am Calibration 
(AmCal, 200 nCi) and active shield trigger signals 
and for absolute timing using 1 pulse per second 
(PPS) pulses from the GPS unit on the gondola. 

We obtain CZT crystals from Redlen Technolo- 
gies 2 , which produce high quality CZT at low cost 
through their patented traveling heater method. We 
bias 64 CZTs at -600V using two EMCO 3 High 
Voltage Power Supplies (HV PS, C12N) mounted on 
the FCB. We apply the HV lead on the cathode us- 
ing a thin (50 /im) Al tape, which is suitable for pres- 
surized balloon flights and convenient for replacing 
DCU units if needed. The Redlen CZTs show < 0.5 
nA/pixel leakage currents under —600V. Each DCU 
is surrounded by Kapton-coated Copper RFI shields 
to minimize the observed interference between the 
units (see [6]). The shields kept the electronics noise 
down to the level of a single DCU operation. 

CZTs were bonded on an interposer board (IPB) 
with two low-cost, low-temperature (< 100 °C) tech- 
niques - initially with low temperature solder bonds 
from Delphon/QuikPak 4 and then with Transient 
Liquid Phase Sintering (TLPS) bonds from Creative 



Table 1 

Telescope Parameters of ProtoEXISTl 



1 http://www.netburner.com/ 

2 http://www.redlen.com 

3 http : //www . emchohighvoltage . com/ 

4 http://quikpak.com 



Parameters 



Values 



Sensitivity (5 a) 
Energy Range 
Energy Resolution 
Field of View 

Angular Resolution 



~ 140 mCrab/hr a 

30 - 300 keV b 

3-5 keV 

19° x 19° (50% Coding) 

9° x 9° (Fully Coded) 

20.3' 



CZT Detector 

Area 

Pixel, Thickness 



64 crystals x (2 x 2 x 0.5 cm 3 ) 
16 x 16 cm 2 
2.5 mm, 5 mm 



Tungsten Mask 12 layers x 0.3 mm thick 

Coding Area 30.1 x 30.1 cm 2 

Pixel, Grid, Thickness 4.7 mm, 0.5 mm, 3.6 mm 

Mask-Det. Separation 90 cm 

CsI(Na) Rear Shield 26 x 26 x 2 cm 3 

Side Passive Shield Semi-graded Pb/Sn/Cu 



241 Am Cal. Source 



200 nCi, 36 cm above the detector 



(a) assuming a 86% QE and a 45% dead time fraction with 
an observed ~ 250 cps background rate in the 30-200 keV 
band (see §4 and Fig. 8b), on-axis and no atmospheric 
absorption, (b) for imaging with >80% mask efficiency as 
limited by mask thickness; the detector energy range is ~ 
30 - 600 keV. 

Electron Inc 5 for final DCUs. TLPS bonding was 
explored for possible future applications and one of 
its main advantages is that the resulting bonding 
is re- workable: the CZT can be detached and re- 
attached if needed. 

The RadNet ASIC was chosen for ProtoEXISTl 
mainly because of its relatively low power consump- 
tion (100 - 150 /iW/pixel) and the expected advance 
of the subsequent ASIC members in the same devel- 
opment family [7,8]. The low power consumption is 
an essential requirement for a large area CZT detec- 
tor with fine pixels, which results in a large number 
of pixels or channels (e.g. 256k pixels in MIRAX). 
In addition, the RadNet ASIC provides a wide dy- 
namic range with low electronics noise suitable for 
the hard X-ray band and allows flexible readout 
modes - multi-pixel pulse profile readout - for depth 
sensing (based on neighbor pixel charges, [9,10]) and 
event reconstruction of multi-pixel triggers arising 
from split charges or Compton scattering [6] . 



http : //creativeelectron . com 



(a) Detector Crystal Unit: 
DCU, 4 cm 2 



(b) Detector Crystal Array: 
DCA.32 cm 2 



(c) Detector Module: 
DM, 256 cm 2 
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Fig. 2. The CZT detector assembly for the ProtoEXISTl imaging detector, (a) Detector Crystal Unit (DCU) consists of 2 x 
2 cm 2 x 0.5 cm CZT bonded via an interposer board to the RadNET ASIC with 1-D 64 inputs (pink), (b) Detector Crystal 
Array (DCA) consists of a 2x4 array of DCUs. (c) Final view of flight Detector Module (DM) covers 16 x 16 cm 2 with a 2 x 
4 array of DCAs. HV bias (— 600V) lines connect to the 4 DCAs on each side through Al tape "jumpers" between DCUs. 



For ProtoEXISTl, we have implemented four 
readout modes for debugging, trigger, normal oper- 
ation, or calibration (see Table 1 in [6]). During the 
flight, we operate the detector in the normal mode, 
where we read signals from eight neighbor pixels 
and two reference pixels in addition to the trigger 
pixels. The total number of neighbor pixels to read 
out varies, depending on the location of the trigger 
pixels: five for edge pixels and three for corner pixels. 
In ProtoEXISTl, we have not recorded the signals 
from adjacent crystals since the electrostatic shields 
prevent the charge induction across the crystals. 

2.2. Mask, Shield and Onboard Calibration Source 

The mask consists of 12 identical stacked layers 
of 0.3 mm thick Tungsten sheets, each of which has 
an identical pattern of a 2 x 2 uniformaly redun- 
dant array (URA) chemically etched out in a 64 x 
64 square pixel configuration. In theory, URAs allow 
a perfect imaging response for a point source with- 
out coding noise. This particular URA pattern used 
for ProtoEXISTl is generated by perfect binary ar- 
rays [12], which allow a square mask pattern ideal 
for our detector geometry (see Fig. 4). The mask 
pixel pitch is a 4.7 mm with 0.5 mm grid to sup- 



port isolated opaque elements. Therefore, the size of 
open elements is 4.2 x 4.2 mm 2 and the total open 
fraction is 40% for normal incidence (the open hole 
fraction is 48%). The mask, mounted on a support 
frame, is located at 90 cm above the detector plane, 
and the URA pattern covers 30 x 30 cm 2 (each mask 
sheet spans 36 x 36 cm 2 for mounting), allowing a 
10° x 10° fully coded field of view (FoV) or 20° x 
20° for 50% coding fraction. 

Graded passive shields made of Pb/Sn/Cu were 
used to limit the X-rays incident through the side of 
the telescope. These side shields are divided into two 
sections (Fig. 1). The upper section sits on top out- 
side of the pressure vessel (PV), mounted all around 
the mask support tower above the top cover of the 
PV to the mask. The lower section is mounted onto 
a similar support tower inside the PV and covers all 
four sides of the detector plane above the rear active 
shield to the bottom side of the top cover of the PV. 
A part of the side shields was not completely cov- 
ered with Sn/Cu, and so the lead sheets were partly 
directly exposed to the detector plane so that fluo- 
rescent X-ray K a and Kp lines from the lead were 
observed in the spectrum (see §4.3). The thickness 
of the side shields are 2.6/0.7/0.4 mm (Pb/Sn/Cu) 
for the lower section, and 2.6/0.3/0.3 mm for the 



Table 2 

Science Data Size and Rate 



Parameters Values 



Event header size a + checksum 48 Bytes 

16 samplings per pixel 32 Bytes 

Raw event size per TV pixel readout 48 + 327V 

Compressed event size per TV pixel readout 48 + 2N 
In normal readout mode, 

Average number of pixels for readout 9.12 pixels 

Average size of raw event data 339.8 Bytes 

Average size of compressed event data 66.2 Bytes 

Measured total event count rate ~ 750 cps b 

Raw data rate 2.2 Mbps 



Compressed data rate 530 kbps 



(a) This does not include the packet header added during 
the transfer between the FCB and the flight computer. 

(b) At 39 km float altitude, and dominated by cosmic-ray 
induced events before shied rejection. X-ray background 
events (30-200 keV) are ~ 250 cps; see §4.2. 

upper section. 

In order to stop the albedo atmospheric X- 
rays and tag charged particle induced background 
events, an active shield of a 2 cm thick 26 x 26 
cm 2 CsI(Na) scintillator crystal from Scintitech 6 
is mounted about 10 cm below the detector plane. 
The Csl scintillator is larger than the active area of 
the detector plane, extending out to the side shields 
in order to reduce the direct open path to the de- 
tector plane from outside. However, the current 
configuration leaves about 4 cm gap between the 
rear and side shields, allowing a high background 
rate during the flight (see §4.2). The optical scintil- 
lation photons are collected by two photomultiplier 
tubes (PMTs, Hamamatsu 7 R1924A) on two op- 
posite sides of the Csl crystal through edge-readout 
wave shifters. Preamps and subsequent electronics 
in a VME crate provide a pulse for tagging shield 
events. A PC- 104 based computer (separate from 
the flight computer below) controls the VME elec- 
tronics. The FPGA on the FCB monitors the two 
discriminator signals from the active shields, and if 
they occur within 2 /isec of the X-ray triggers, the 
X-ray event is tagged. 

In order to monitor the gain variations of the de- 
tector, a 241 Am calibration source (AmCal) was as- 
sembled and mounted. The AmCal source consists of 



a 200 nCi 241 Am doped in a plastic scintillator (Isop- 
tope Products Laboratories 8 ), which is coupled to a 
PMT (Hamamatsu R7400U). The scintillation light 
induced from the a particle coincident with a 60 keV 
X-ray is collected by the PMT and amplified through 
a separate channel in the shield VME electronics for 
tagging the Am events. In addition to the onboard 
calibration source, we can command all of the 64 
RadNET ASICs to inject on-chip pulses into every 
pixel at a configurable frequency, and we read them 
out in the same way as normal X-ray events. This 
monitors the dead time and the electronics gain vari- 
ation of the readout system throughout the flight. 

The fully assembled detector plane (Fig. lc) is 
mounted on an Al frame and surrounded by a thin 
copper box on the four sides and the rear side, which 
forms a detector shield enclosure. The top cover of 
the enclosure box is a 150 jiva thick Al sheet. Two 
fans are mounted on the side of the enclosure and 
another two on the bottom for temperature control 
through air (N2) circulation through the PV. The 
fan power is electrically isolated from the power for 
the detector plane to minimize the electronics noise. 

2.3. Flight Computer, Data Rate and Transfer 

The flight computer is a 2 GHz Intel dual core 
compact PC (AEC-6920 by AAEON Technology 9 ) 
with a 128 GB solid state disk designed to operate in 
extreme environments. The flight computer collects 
the science and HK data from the FCB through Eth- 
ernet links, stores them on the onboard hard disk, 
and sends a compressed version of the data to the 
transmitter (see below). It also issues to the FCB 
the commands uploaded through a RS422 port from 
the gondola computer [5]. In order to handle these 
tasks, the flight computer runs three basic programs 
- a commmand listener (CL) program, a data ac- 
quistion (DAQ) program, and a transmitter handler 
(TH) program under a linux (Fedora 8) operating 
system. A cron job based watchdog script monitors 
the health of these three programs, and terminates 
or relaunches them automatically if needed. 

Table 2 summarizes the data size and measured 
rate (see also Fig. 6). The raw data rate is about 2200 
kbps for 750 cps. The DAQ program saves the en- 
tire raw data stream on the onboard solid state hard 
disk. In order to meet the telemetry limit (400 kbps), 



6 http://www.scintitech.com 

7 http://www.hamamastu.com 



8 http://www. isotopeproducts.com/ 

9 http://www.aaeon.com 
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the DAQ program also compresses the data by calcu- 
lating a single pulse height value from 16 samples of 
the pulse profile (see Fig. 2 in [6]), and sends only the 
resulting pulse height instead of the full 16 samples 
to the TH program. The corresponding compressed 
data rate is about 530 kbps, which is about 24% of 
the raw data rate instead of 8% (l/16th) due to the 
overhead of the event header. The compressed data 
rate can still exceed the telemetry rate. Therefore, 
the DAQ program monitors the data rate, and drops 
a fraction of events randomly (only for telemetry) to 
meet the telemetry bandwidth so that the full data 
stream for a given event can be transmitted to the 
ground during the flight. The data reduction frac- 
tion for telemetry is continuously adjusted during 
the flight to keep the compressed data rate below 
300 kbps, which allows a 100 kbps margin for the 
packet header, the HK data or possible variations 
of the telemetry limit. For debugging purpose, the 
DAQ program can also send the raw data without 
compression at request. Transmission modes (raw 
or compressed) and the data reduction fraction for 
the telemetry can be controlled by the commands 
uploaded from the ground. 

In addition to the science data, the HK data and 
the command echo are merged and transmitted to 
the ground without any compression since their rate 
is relatively small. There are four major HK data 
streams: the FCB HK data reporting the various 
input voltages and currents including the HV PS 
and the temperature of the FCB, the shield HK data 
reporting the active shield and AmCal count rate 
and status of the HV PS for PMTs, the DAQ HK 
data reporting the status of the flight computer such 
as free hard disk space, and the Power HK data 
reporting the power input voltages and the currents 
coming into the PV and temperatures and pressures 



of a few points inside the PV. On the ground, we 
monitor and save the transmitted data, and issue 
commmands accordingly. 



3. Pay load Integration and Test 

The detailed description of the payload integra- 
tion is found in [4,5]. Here we briefly summarize the 
integration and test of the detector plane. 

The CZT detector plane for the ProtoEXISTl 
telescope was fully assembled (Fig. lc) at Harvard 
University on 2009 August 31. In September, 2009, 
the detector plane was mounted inside the PV along 
with the side passive shields, the Csl rear active 
shield, the AmCal source and the electronics. Then, 
the ProtoEXISTl telescope payload - the PV and 
the mask tower - was integrated into the Harvard 
gondola now including a new pointing and aspect 
control system consisting of a daytime star camera, 
differential global positioning system (DGPS) units, 
and a magnetometer (Fig. 3) [4,5]. 

After the payload integration, we irradiated the 
X-ray telescope using a ~ 5 mCi 57 Co radioactive 
source. The source was positioned about 6.5 m away 
from the mask within a few degrees of the X-ray tele- 
scope axis. Fig. 4a shows the detector count map of 
the source in the 80 - 150 keV band from about 15 
min of the irradiation of the 57 Co source. For illus- 
tration, the 2x2 URA mask pattern is shown in 
Fig. 4b with a section highlighting the matched pat- 
tern observed in the detector image. Fig. 4c shows 
the reconstructed sky image of the fully coded FoV 
using a balanced cross-correlation (e.g. see [11]) be- 
tween the mask pattern and the detector image and 
accounting for magnification of the shadow pattern 
due to the finite source distance. The measured S/N 
(~ 80) of the source is a bit lower than the Poisson- 
noise based S/N because of the coding noise intro- 
duced by incomplete sampling (91%) of a full cycle 
of the URA pattern due to magnification. 

For energy calibration, we use the pre-flight radi- 
ation data using the 57 Co source. For a given pixel 
at a given trigger capacitor, we identify the raw 
ADC values of three known energies - zero keV from 
the reference pixel measurements, 122 keV from the 
57 Co source and 650 keV from the internal pulser 
inputs - by a fit using a gaussian function to the raw 
spectra. The conversion from other raw ADC values 
to energies uses a 2nd order polynominal fit between 
the raw ADC values and X-ray energies of the three 
known values [9] . After this relatively quick calibra- 



(a) 




(b) 




Fig. 4. (a) The detector count map generated by a 15 min irradiation of an external 57 Co source at 6.5 m away from the 
telescope, (b) The 2x2 URA mask pattern: the (orange) shaded section of the mask pattern indiciates the coding area of this 
particular configuration, (c) The reconstructed sky image in the fully coded FoV using 5.35 mm magnified mask pixel pitch. 



tion, the detector shows about 4 keV FWHM at 60 
KeV lines and 10 keV at 120 keV. The latter indi- 
cates that there is room for improvement in the cal- 
ibration (see [9,6,5]): e.g. the pulser is a less reliable 
energy calibrator, and the charge split or depth cor- 
rection is not performed. 

The tagging efficiency of the AmCal events in the 
CZT detector is measured about 74% (see §4.2), and 
does not increase with the coincidence window as 
long as the window is longer than 2 /isec. Therefore, 
the other 26% AmCal events are suspected to fail to 
generate a tagging signal all together probably due 
to insufficient light collection in the PMT from the 
5 MeV a particle of 241 Am decays near the edge of 
the scintillator. The tagging efficiency of high en- 
ergy charged particles by the Csl shield is estimated 
to be around 50% on the ground, based on the high 
energy events induced by atmospheric muons, which 
are recorded near the upper limit (^ 1 MeV) of en- 
ergy deposit for a pixel. This is roughly consistent 
with the solid angle coverage of the rear shield. This 
tagging efficiency is likely energy dependent. 

4. Flight Performance 

4.1. Flight Overview 

After 6 previous launch attempts, early in the 
morning on Oct 9, 2009, all the systems were 
checked twice with an external power supply for 
flight, and then before balloon inflation began, they 
were switched to the onboard battery, waiting for 
launch. The procedures for the initialization se- 
quence of the detector are, in order: powering up the 
system, activating the internal ASIC pulsers, dis- 
abling the known hot pixels, applying the HV bias 
on the CZT crystals at —600 V (gradually ramped 



up over ~ 10 min), and setting the DCU thresholds 
gradually down to ~ 30 keV over ~ 15 min. Each 
DCU has its own unique threshold setting but all 
64 pixels in a DCU share a common threshold. The 
complete sequence takes about 30 min and each 
procedure can be either manually or autonomously 
executed. The internal ASIC pulsers were run in all 
of the 64 DCUs at 0.3 Hz throughout the flight for 
a dead time calibration. The active shield system is 
also initialized right after power- up. 

The pay load was launched at 14:30 UT by a 40 
million m 3 high altitude balloon from the NASA 
CSBF at Ft. Sumner, New Mexico. Fig. 5 shows (a) 
the inflated balloon bubble (top ~ 1/4 of the bal- 
loon) right before the launch, (b) during the liftoff 
and (c) the flight path on a Google map. Fig. 6 shows 
the aspect and pointing history of the system dur- 
ing the flight and the raw count rate history of sev- 
eral key event types recorded in the detector system. 
The payload reached the float altitude of 40 km at 
17:00 UT (blue in Fig. 6a). The flight duration was 
only about 7.5 hr with 5.5 hr at float altitude due to 
strong winds (~60 knots) and it was terminated near 
Hayes, Kansas, due to both the ~ 500 km teleme- 
try limit and the NASA safety regulations that pro- 
hibit flights over densely populated area. The land- 
ing was relatively harsh due to high surface winds, 
and the gondola frame sustained some damage, but 
the payload was recovered without any damage. 

When the payload reached the float altitude, the 
Crab pulsar, one of a few sources bright enough for 
detection in a short observation for a hard X-ray 
telescope of this moderate size, had just set. Until 
Cyg X-l rose about 4 hr later, we succesfully per- 
formed a number of system checks such as shield 
threshold tests (18:00 UT), while attempting track- 
ing 3C 273. From 21:00 UT, we observed Cyg X- 




Fig. 5. The inflated balloon "bubble" (top ~ 1/4 of balloon) right before the launch (a), the liftoff (b) and the flight path on 
a Google map (from Ft. Sumner, New Mexico to Hayes, Kansas) (c). The inset in (b) shows the payload during the liftoff. 



1 until the termination of the flight. After an hour 
observation of Cyg X-l, we reset the power of the 
detector system at 22:00 UT for a test. The system 
came back up without a glitch. After the power re- 
set and the initialization sequence, we had less than 
a 10 min of fine-pointed observation of Cyg X-l be- 
fore we were requested to terminate the flight. 

The X-ray telescope performed excellently 
throughout the flight, but the aspect and pointing 
system suffered a few major issues: until in the 
middle of the Cyg X-l observation, the telescope 
could not track a target properly. First, we lost the 
DGPS units, possibly due to insufficient thermal 
insulation, at around 15:46 UT during the ascent. 
Second, a malfunction of the star camera software 
disabled the live feedback to the pointing system 
and forced us to rely on the magnetometer and 
the elevation axis for target tracking. The absolute 
error of the raw azimuth values reported by the 
magnetometer can be as large as 10° without cal- 
ibration, depending on the orientation. Third, the 
elevation axis failed to follow the target properly 
until stablized later in the flight (21:55 UT), likely 
due to insufficient torque for the motors. 

Despite these issues, Cyg X-l was within the par- 
tially coded FoV above 60% coding fraction and it 
was mostly < 2-3° off in the elevation axis and 6- 
7° off in azimuth from the X-ray axis (see Fig. 12). 
We calibrated the raw azimuth values using those 
reported by the DGPS units during the accent be- 
fore they failed when the payload was freely spin- 
ning (14:30 - 15:46 UT; black line in Fig 6a). This 
calibration reduces the absolute error of the azimuth 
down to ~ 1-2°. Fortunately, a set of the sky im- 
ages captured by the star camera were later identi- 
fied and used to verify and correct the aspect and 



pointing history derived by the magnetometer and 
the elevation axis. The time stamps of these images 
are marked by (green) 'X's at the right-most side in 
the top panel of Fig. 6. 

4.2. Event Rate 

Fig. 6b shows the raw count rate history of the 
detector system during the flight. The total event 
rate (black) is the total trigger rate in the detec- 
tor system, which is dominated by the cosmic ray 
induced interactions. The shield tagged event rate 
is shown in red, and the X-ray event rate (blue) is 
from the events without a shield tag. The AmCal 
rate (purple) is the event rate with an AmCal tag. 
The pulser rate (green) is the internal ASIC pulser 
rate. Note that each raw event is recorded in a unit 
of a DCU pair, counting a multi-pixel trigger event 
as one event but counting a multi-DCU-pair trigger 
event as multiple events. The latter means that the 
0.3 Hz pulser injected over the entire detector plane 
should be recorded with rate 0.0394 counts s _1 cm -2 
when there is no dead time 10 . 

The absolute timing accuracy of events is in an 
order of 0.5 /is, which is determined by the sampling 
frequency of the RadNET ASIC. We save the time 
tag of each event with a 0.1 ms resolution. A typical 
readout time for an event with a 10 pixel readout 
is about 5 ms, which leads to a dead time for event 
readout in a given DCU pair, but not for the rest of 
the DCU pairs since each DCU pair operates inde- 
pendently. The analysis show that almost all of the 
X-ray events with multi-pixel triggers in a DCU pair 



10 0.0394 = 0.3 Hz x 32 DCU pairs / 243.4 cm 2 , where the 
detector area is 1.95 x 1.95 cm 2 x 64 DCUs. 




Cal events due to missing tags and is not corrected 
for dead time. To extract the proper X-ray event 
rate and spectra, we use the following relations. 
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Fig. 6. (a) The aspect and pointing history of the flight. The 
azimuth angles (black) are based on the DGPS-calibrated 
magnetometer, and the elevation angles (red) are based on 
the elevation axis. The altitude (blue) is recorded by external 
pressure gauges. The (green) 'X's mark the azimuth and 
elevation values drived by a set of sky images captured by 
the star camera. The top axis labels the key events of the 
flight, (b) The raw detector count rate history of the flight. 
From top, the total event rate (solid black), the X-ray event 
rate (blue), the shield tagged event rate (red), the AmCal 
tagged event rate (solid blue), and the pulser event rate 
(green). The dip in pulser rate at 15-16 hr is due to increased 
dead time at the high count rates encountered in the Pfotzer 
maximum in balloon ascent. 

are truly simultaneous (not a random coincidence 
of multiple events), whereas about 9% or less of the 
X-ray events with multi-DCU triggers are truly si- 
multaneous. In the case of multi-pixel triggers, al- 
most all of the X-ray events are four pixel triggers 
or less. Therefore, for imaging and spectral analysis 
of Cyg X-l and the AmCal events, we mainly focus 
on the low energy X-ray events with four pixel trig- 
gers or less, and we consider all the low energy X- 
ray events with multi-DCU- pair triggers as indepen- 
dent events. For the charged particle induced events 
with a shield tag, about 30% triggered multiple DCU 
pairs. Many of them are truly simultaneous, gener- 
ating large tracks (Fig. 7), which illustrates a poten- 
tial of large area CZT imaging detectors as cosmic 
ray detectors. 

The raw X-ray event rate in Fig. 6b contain Am- 



Ot = fuve(RAm + Rx) ~ ^Am 
^Am = /Live(/AT^Am + /rT^x) 



(i) 

(2) 



where /Live? /rt, /at are the live time fraction, the 
random coincidence fraction by AmCal tag and the 
AmCal tagging efficiency respectively, and rx and 
rAm are, respectively, the measured raw count rates 
of the X-ray events (no shield or AmCal tag) and the 
AmCal events (an AmCal tag, but no shield tag); 
Rx and Raih are the corresponding true count rates 
(i.e. the incident rate if the detection efficiency is 
100%). Fig. 8 shows these corrected rates (Rx and 
RAm) for various cases and Fig. 9 for the energy 
spectra before and after correction. 

The live time fraction (/Live) is estimated based 
on the captured internal pulser fraction and ran- 
dom coincidence rate of AmCal tagging (/rt) is es- 
timated based on hard X-ray events in the 80-200 
keV band where all the AmCal tagged events are 
due to random coincidences. (Fig. 8a). Pulser events 
are easy to identify due to simultanoues injection of 
650 keV equivalent signals into all 128 pixels of each 
DCU pair, which produce a narrow line around 83 
MeV (for summed pixels) in the pulse height his- 
togram, (e.g. green in Fig. 9c). The pulser recov- 
ery fraction (= /Live) is about 97% on the ground, 
and it dropped down to 40% during the ascent due 
to a large increase of background, and then stayed 
around 55-60% during the flight. n Fig. 8a (red) 
shows /rt stayed fairly constant throughout the 
flight. For the rest of the calculation, we use the av- 
erage value of /rt, 3.6%. 

We set /at so that the resulting energy histogram 
of the X-ray events free of the AmCal events shows 
no positive or negative bump around the 60 keV 
(black curve in Fig. 9b). We did this for both the 
flight spectra and the prelaunch ground spectra, and 
for /Live, we use the average value of each data set. 
In order to have a smooth continuum at 60 keV, the 



11 The event readout time in a DCU pair ranges from ~ 5 to 
25 ms, depending on the number of readout pixels. For an 
event occurance rate of 1500 cps over the full detector plane, 
which is equivalent to about 47 cps per a DCU pair, the 
event capture fraction under an on-average ~ 15 ms event 
readout time would be about 50% (= exp(-47x 1.5 x 10 -2 )). 
Therefore, the observed total event rate of ~ 750 cps is 
consistent with the ~ 50% live time fraction. 
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Fig. 7. Examples of multiple DCU trigger events induced by cosmic ray interactions during the flight. The number of the 
simultaneously triggered pixels and the total energy deposited by the events are labeled. The color scale shows the energy 
deposit in each pixel in keV. Since the energy deposit is saturated at ~ 1 MeV, the recorded total energy deposit is a lower limit. 



allowed range for /at is between 0.735 and 0.740; 
we therefore used 0.738 for /at for the calculation. 
Fig. 8b shows the corrected total event (gray) , X- 
ray event (black) and AmCal (red) event rates. The 
total event rate (grey) reached about 9 counts s _1 
during the ascent and dropped to 4.5 counts 



cm 

s _1 cm -2 at the float altitude and gradually in- 
creased to 5.4 counts s _1 cm -2 later in the flight. 
The X-ray event rate (black) stayed about 2.3-2.7 
counts s _1 cm -2 at the float altitude. Since the trig- 
ger rate in the Csl shield (~ 100 keV threshold), 
which likely represents the charged particle induced 
event rate, (not shown in the figure) was 2.1 - 2.7 
counts s _1 cm -2 , the combined rate of the X-ray and 
charged particle induced events makes up the ob- 
served total count rate of 5.4 counts s _1 cm -2 . How- 
ever, the X-ray event rate still contains a large frac- 
tion of charged particle induced events, considering 
the relatively low shield tagging efficiency (~50%, 
§3), a higher threshold of the Csl shield, and sec- 
ondary background due to interactions in the sur- 
rounding side shield and mounting structures. To 
resolve the precise composition of the X-ray back- 
ground rate, detailed Monte-Carlo simulations of X- 
ray and charged particle transport in the telescope 
geometry are needed. 

The low energy X-ray events (30-200 keV, blue) 
stayed at around 0.86 -1.1 counts s _1 cm -2 . If we as- 
sume a shield tagging efficiency of 50% (§3), the low 
energy X-ray event rate was about 0.5 - 0.7 counts 
s _1 cm -2 (purple). Assuming an albedo X-ray flux 
of 13.MT 1 - 81 ph cm" 2 s" 1 keV" 1 at 3.5 g cm" 2 [13] 
and an aperture FoV of 0.1 sr (FWHM) with a 40% 
open fraction, we expect the aperture X-ray back- 
ground rate of ~ 0.04 counts s _1 cm -2 in the same 
band. In addition, the gap between the rear and side 



shields allows a projection-corrected FoV of 0.5 sr, 
which can contribute ~ 0.45 counts s _1 cm -2 . We 
also estimate a total flux of ~ 0.1 counts s _1 cm -2 
in the Pb fluorescent X-rays. The combined aper- 
ture and leakage X-ray flux in the 30-200 keV band 
is estimated to be ~ 0.6 counts s _1 cm -2 , which is 
consistent with the observed X-ray event rate. This 
estimate lacks some details such as attenuation by 
surrounding structures for X-rays coming below, the 
contribution from the internal background or an in- 
crease of the active area due to gaps between CZTs 
but their effects are rather small or their contribu- 
tions cancel each other to some extent. 

The corrected AmCal rate (i?Am> red) stayed rela- 
tively constant at 0.128 ± 0.006 counts s _1 cm -2 , in- 
dicating our analysis works well. Based on the known 
source strength of the AmCal sources with ~ 2% at- 
tenuation from the detector enclosure window, etc, 
one can calculate the overall detection efficiency at 
60 keV to be 0.86 ± 0.04; the - 10% loss of QE is 
due to a variety of effects (dead surface layer, etc) 
and will be reported for several energies in [5]. The 
photo-peak efficiency in the 55 - 65 keV range is 0.54 
± 0.02. The error values are only of statistical origin, 
and the simplified geometric assumption of a point 
source can add an additional 5% systematic error. 

4.3. Detector Spectra 

Fig. 9 shows the energy histograms of various 
event types before and after correction by Eq. 1 and 
2 in §4.2. The top panels show the low energy sec- 
tion of the histogram in photon flux, and the bot- 
tom panels show the full range of the histogram in 
energy flux unit. Fig. 9d breaks up the X-ray events 
(no shield tag, no AmCal tag) by the number of 
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Fig. 8. (a) The pulser capture fraction (= the live time frac- 
tion, black) and random tagging fraction of the AmCal source 
events (red), (b) The count rate history of the corrected data 
by Eq. (1) & (2). Note the correction fixes only the untagged 
AmCal event rate, so that the X-ray event rate still includes 
contributions from the charged particle induced events. 



triggering pixels. The flight data in the figure are 
from the events recorded between 17:50 and 21:53 
UT, and the ground data between 13:35 and 14:40 
UT. In the case of the corrected AmCal spectra the 
ground data (not shown) agree with the flight data 
(red in Fig. 9b & c). Note that, after correction, the 
AmCal spectra lack high energy events (>65 keV). 
A 5 mm thick CZT detector is about 80% trans- 
parent to 500-600 keV photons, and the higher en- 
ergy events in the spectra are mostly from the multi- 
pixel triggering events such as the charged parti- 
cle induced events and the pulser events at 83 MeV 
(green) (Fig. 9c). The peaks at 75 and 85 keV are 
due to the Ka and K/3 lines of the lead fluorencent 
X-rays as aforementioned. A mild bump around 511 
keV originating from electron positron pair produc- 
tion is visible in the single pixel trigger event (blue) 
of Fig. 9d. The origin of the peak at ^ 1 MeV is not 
cosmic but is due to the upper limit of the dynamic 
range of the ASIC preamp channel (i.e. the maxi- 
mally allowed energy deposit for a pixel), and the 
multi-pixel readout generates the additional peaks 
at its multiples around 2 and 3 MeV. 



4.4. Temperature Dependent Gain Variation 

The detector system is mounted inside of the PV, 
which is pressurized with dry N2 and thermally con- 
trolled by a set of heaters and fans. The detector 
module experienced a ~ 10° C temperature variation 
during the flight. Fig. 11 shows the AmCal spectra 
summed over the entire detector plane as a function 
of time in comparison with the temperature vari- 
ation (blue). For easy comparison, the energy his- 
togram is renormalized at each time bin. The figure 
shows that the detector gain correlates with the tem- 
perature (AE (keV) « 0.116AT (°C)). The tem- 
perature dependent gain change is likely due to the 
capacitance variation of sampling capacitors in the 
ASIC, and the CZT properties do not change ob- 
servably over this temperature variation. A similar 
trend is observable with the pulser run. Unlike the 
gain, the energy resolution does not vary much with 
the temperature. For more information on the tem- 
perature dependent detector response of the Pro- 
toEXISTl detector, see [5]. 

5. Imaging Cyg X-l 

We observed Cyg X-l for about 100 mins at the 
end of the flight. After an hour into the observa- 
tion, which had large pointing excursions (mostly in 
elevation; see Fig. 12a), we performed a power cy- 
cle test of the detector system, which took about 30 
min. Thus we have a total of an hour and ten min- 
utes of usable data with only the final 10 min having 
relatively stable (± 10') pointing (see Fig. 12). Here 
we decribe the aspect calibration, image processing 
and detection of Cyg X-l as well as its measured 
energy spectrum. 

5.1. Image Reconstruction 

The event distributions in the detector observed 
during the flight are substantially different from 
those measured on the ground. For instance, the 
dominance of charged particle interactions increases 
the fraction of multi-pixel events, which more likely 
produce an artificial pattern of event distribution. 
In low energy X-ray events with < 4 pixel triggers, 
there are a few factors causing the non- uniformity 
in the event distribution in addition to the usual 
pixel or crystal dependent efficiency variation. 
Fig. 10a shows an example raw detector image of 
single pixel trigger events. In order to visualize the 
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Fig. 9. (a) The energy-calibrated observed flight spectra classified by shield and AmCal tags, (b) & (c) the corrected flight 
spectra by Eq. 1 &; 2 in comparison with the corrected ground X-ray spectra (blue), and (d) the corrected flight spectrum of 
the X-ray events and their composition by the (labeled) number of triggering pixels. The left panels show the spectra in the 
photon flux and the right in the energy flux with a wider energy range. In (a), the X-ray spectra (rx) represent events without 
an AmCal or shield tag, but they still contain the AmCal events, as indicated in a large peak at 60 keV. After correction in 
(b) and (c), the X-ray spectra (Rx) do not have any contribution from the AmCal source. In (c) the X-ray spectra of the 
corrected flight data (black) are substantially brighter than those of the ground data (blue), but the AmCal and pulser spectra 
of the corrected flight data are nearly identical to those of the ground data (not shown). The vertical lines in (d) mark 511 
keV and the maximally allowed energy deposit in a pixel (~ 1 MeV); higher energy peaks are due to 2 or 3 pixel pileup. 



non- uniformity more clearly, we also show the event 
distribution folded onto a DCU pair (i.e. all 32 DCU 
pairs are stacked), as illustrated in the small panel. 
In this example, the first noticeable non-uniformity 
is the count enhancement of the edge pixels due to 
gaps (900 /im for 1.95 cm x 1.95 cm crystals) be- 
tween CZT crystals that provide additional surface 
area for accepting off- axis X-ray events. In addition, 
the current sequential readout scheme performed 
in an increasing order of pixel number makes the 
data of the last few pixels more susceptible to 
corruption. In summary, non-uniformities in the 
detector image have various origins: for some, the 
cause is obvious and for others, not. For imaging, 
we rely on background subtraction to eliminate the 
non-uniformities regardless of their cause. We also 
explored an option of using a synthetic background 
pattern for subtraction, which was generated based 
on the observed non-uniformities. The synthetic 



background pattern does not work as well, which is 
likely due to unknown and thus missing components 
of the non-uniformities in the synthetic pattern. 

We generated a background detector plane im- 
age using the data set between 18:30 and 20:58 UT, 
which is just before pointing on Cyg X-l and so 
the closest in-flight background free of bright X-ray 
sources. Then we subtracted the scaled (for total 
counts) background image from the detector plane 
images of the observation. Note the overall back- 
ground rate varies with altitude and elevation, but 
the pattern of non-uniformities appears less sensi- 
tive to altitude or elevation as demonstrated by the 
uniformity of the background subtracted images be- 
low (Fig. 10b). Since the background dominates the 
source counts (< 5%), the scale factor used for back- 
ground subtraction is simply the ratio of the total 
count of the observation to the total count of the 
background image, which makes the total counts of 
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Fig. 10. Example detector plane images for single pixel trigger events in 30 - 200 keV: (a) before background subtraction, the 
count varies from pixel to pixel more than 10cr above the average and (b) after background subtraction, the variation is within 
~ 3cr. The DCU-pair folded event distribution in the left top of each panel reveals the non-uniformity more clearly before the 
background subtraction: e.g. an excess of edge pixel counts and a deficit in a pattern that resembles the number '11'. 



the subtracted image zero. The detector had 37 dead 
or disabled pixels (out of 4096 total), which include 
two pixels disabled during the flight. Even after dis- 
abling these hot pixels, there are a few pixels (usu- 
ally less than five) showing about a factor of 5 or 
10 larger count rate than the rest (8- 10a above the 
average). In a detector plane image, we set five 12 
highest pixel counts to the next-highest value. 

Fig. 10b shows an example of the background 
subtracted detector plane images using single pixel 
trigger events in the combined set from the Cyg X- 
1 observation. The background subtraction elimi- 
nates most of the outstanding non-uniformities. For 
the Cyg X-l observation, we generate such a back- 
ground subtracted detector image every 12 s from 
20:57:38 to 21:53:02 UT before the system reset (the 
power-cycle test) and from 22:31:08 UT to 22:39:32 
after the reset. We produce the corresponding re- 
constructed sky images using a balanced correlation 
[ref] with the mask pattern using FFT. The detec- 
tor image of 64 x 64 pixels is rebinned into 8192 
x 8192 pixels in order to properly scale the mis- 
match of the mask (4.7 mm) and detector pixel pitch 
(5.1/2 mm) 13 . Each sky image (every 12 s) is then 
rebinned to 768 x 768 pixels covering down to the 
10% coding fraction, aspected corrected and then 



12 The final results are not sensitive to the exact choice of 
this number (from 5 to ~ 10) as long as the extreme outliers 
(8-10cr above the average) are removed. 

13 The pixel pitch in a CZT crystal is 2.46 mm, but the gaps 
between the detectors make the overall pixel pitch equivalent 
to 2.55 mm (the detector pitch is 20.4 mm). 



summed to produce the final image, weighted by the 
coding fraction. 

5.2. Aspect Correction 

Fig. 12a shows the aspect and pointing history 
of the Cyg X-l observation. The star camera im- 
ages are solved using astrometry.net [14] after post- 
processing the raw images to remove several streaks 
due to a defect in the CCD (see also [4]). The star 
camera images with a good aspect solution were only 
available for a small fraction of the observation (5 
images before the reset around 21:15 UT and 63 im- 
ages after the reset). For the pre-reset data, most of 
which do not have a matching star camera image, we 
correct the X-ray images using the azimuth and el- 
evation given by the DGPS calibrated magnetome- 
ter and the digital readout of the elevation axis with 
shaft angle encoder (SAE) resolution of sub arcmin. 
For about 10 min of the usable post-reset data, we 
use the apect information by the star camera images 
for the relative correction of the X-ray images. 

Fig. 12b shows the relative offset between the as- 
pect of the star camera and the aspect information 
calculated by the magnetometer and elevation axis. 
During the post-reset observation, the relative off- 
set is within 20'- 30' in elevation and < 10' in az- 
imuth, but there seems to be a long term drift in 
the azimuth between the pre and post reset data, 
which can be as large as 50'. This long term drift is 
the residual offset after the DGPS calibration, and 
it originates from the azimuth hysteresis of the mag- 
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Fig. 11. Monitoring the gain variation during the flight using 
AmCal sources. The normallized energy histogram map as a 
function of time overlaid with the temperature of the FCB. 
The gain change AE (keV) « 0.116AT (°C). 



netometer, which can be as large as a degree or two. 
In order to compensate for this offset, the aspect of 
the X-ray image around 21:15 is corrected to match 
the post-reset data, and the rest of the pre-reset data 
is corrected accordingly. 

A caveat of the aspect correction described above 
is that the offset between the X-ray axis and the op- 
tical axis (or the axis of the star camera) must be 
known. A delay in the payload integration left no 
time for the X-ray vs. star camera boresight mea- 
surement before the launch. In order to resolve an 
ambiguity in aspect correction due to the unknown 
boresight, we first correct the pre-reset data accord- 
ing to the aspect information reported by the mag- 
netometer and the elevation axis only relative to the 
image at 21:15, assuming there is no X-ray boresight 
offset. The combined image shows a point source of 
about 6.5a at about 6° off the center of the X-ray im- 
age in the 30-100 keV band. Assuming this is CygX- 
1, one can estimate the X-ray boresight, which was 
calculated to be 31'. This X-ray boresight value is 
reasonable, given the uncertainties in mounting the 
X-ray telescope and the star camera. 

Assuming the 31' X-ray boresight offset, we re- 
apply the aspect correction on the proper X-ray axis, 
and recombine the X-ray images. In the recombined 
image of the pre-reset data, the same point source 
now appeared at 6.8a, indicating the proper bore- 
sighting has increased the signal. Finally we added 
the 10 min of the post-reset data using the aspect 
information given by the star camera images. In the 
final image, the point source appears at 7.2<j. This 
solution gives Cyg X-l at the expected location with 
a reasonable boresight of the X-ray telescope to the 
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Fig. 12. (a) The aspect and pointing history and (b) the 
relative offset of the aspect given by the star camera vs. as- 
pect given by by the magnotometer and the elevation axis 
shaft angle encoder during the observation of Cyg X-l. The 
detector system was reset at around 22:00 UT for a test of 
in-flight recylcing of power and initialization. For imaging 
and spectral analysis of Cyg X-l, the data taken from 20:57 
to 21:53 UT (pre-reset: marked "on Cyg X-l" to "HV off') 
and from 22:31 to 22:39 UT (post reset: marked "Reinitial- 
ized" to "Stowing") were used. 

star camera. Fig. 13a shows the reconstructed sky 
image in the 30-100 keV band. The X-ray and opti- 
cal axis are marked by '+' and V. The offset of Cyg 
X-l from the X-ray axis is 6.4°, resulting in a partial 
coding fraction of 81% and thus 19% vignetting. 



5.3. X-ray Spectrum of Cyg X-l 

In order to calculate the X-ray spectrum of Cyg 
X-l, we repeat the above imaging procedure in the 
30-50, 50-80, 80-125, and 125-190 keV bands. The 
thumbnail images in Fig. 13b show the reconstructed 
sky images around Cyg X-l in these energy bands. 

Due to the relatively large uncertainties (~ 30') 
of the aspect information given by the magnetome- 
ter and the elevation axis, a point source would ap- 
pear somewhat broadened in the reconstructed sky 
images. For simplicity, we assume the source in the 
image is a point source and so the maximum of the 
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Fig. 13. Imaging Cyg X-l: (a) the reconstructed sky image in the 30-100 keV band above 60% coding fraction, (b) the 2° x 
2° thumbnail images around Cyg X-l in the 30-50, 50-80, 80-125 and 125-190 keV bands. In (a), '+' indicates the center of 
the X-ray image, and 'x' indicates the star camera axis. A ~ 4cr blob around (R.A.,DEC.) = (19h 40m, 26°) is a ghost image of 
Cyg X-l, caused by the cyclic nature of the mask pattern (URA). The noise around (19h 36m, 37°) is likely due to incomplete 
aspect correction, in particular, large elevation jitters, since they are aligned along the elevation direction from Cyg X-l. The 
red circles in (b) show the angular resolution (20 ; diameter). 



" blob" in the source image gives the detected source 
counts. Using the maximum source counts is valid in 
our imaging procedure when the background counts 
are much greater than the source counts. In reality, 
this will give us a lower bound for the source counts. 
We calculate the source counts in each band, 
and apply a series of corrections to get the incident 
rate, which include the average dead time (45%), 
the mask open fraction (38.7%), the average cod- 
ing fraction (81%), the average coded mask auto- 
collimation effect (2.6%), the energy-dependent QE 
and atmospheric absorption. We assumed the QE 
of 86% in the 60-100 keV range and for the rest of 
the energies, we use a linear interpolation assum- 
ing 54% at 30 keV (where only photo peaks are 
collected with a ~ 30 keV threshold) and 68% at 
200 keV (based on the total photon cross-section 
of CZT). For atmospheric absorption, we use the 
atmospheric depth of 3.5 g cm -2 at the average 
altitude 38.5 km [15] with the average elevation of 



42.7°. In order to allow energy dependence, we use 
the attenuation of photons in the atmosphere in 
[16] (p 203). For both QE and atmospheric absorp- 
tion, we calculate the energy bin averaged estimates 
using a power law spectrum of photon index T=1.7. 
The resulting spectrum is shown in Fig. 14. The 
error bars along the x and y-axis represent the en- 
ergy bin size and the statistical errors respectively. 
The red point shows the range of the source flux in 
15-50 keV measured by Swift /BAT 14 within 12 
hours of the ProtoEXISTl observation. The closest 
BAT observations are about 1-2 hours before or af- 
ter the ProtoEXISTl observation. The dashed line 
shows a representative hard X-ray spectral model 
of Cyg X-l in the literature [17]. The observed flux 
is a bit lower than the other measurements but it 
is certainly in an acceptable range, given the high 



14 http : //heasarc . gsfc.nasa.gov/docs/swift/results/ 
transients/ 
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ASIC, developed for NuSTAR 15 will be used in the 
8x 8 DCU close-tiled CZT detectors of ProtoEX- 
IST2. The Nu-ASIC is developed in the same ASIC 
family line as the RadNET ASIC, so they share a 
similar basic operational architecture, which simpli- 
fies the transition. We will combine the Nu-ASIC 
front-end readout electronics with a new packaging 
and modularization architecture derived from the 
ProtoEXISTl system. For more details of the Pro- 
toEXIST2 detector development, see also [18]. 

7. Acknowledgement 



Fig. 14. The observed hard X-ray spectrum of Cyg X-l. 
The (red) data point represents the range of the observed 
flux in the 15-50 keV band by the BAT within 12 hours 
of the ProtoEXISTl observation. The broken power model 
(dashed) is from Steinle, et al, 1982 [17]. The solid line is a 
power law fit (~ E~ r ) to the data points, where the photon 
index r = 1.7 ± 0.5. 



X-ray variability of the source. The observed spec- 
trum agrees with the known spectral models within 
3a and it can be fit by a simple power law model 
with photon index Y—l.l ±0.5 (the solid line). The 
reduced x 2 of the fit is small (~ 0.02) due to the 
large uncertainties. The observed spectrum is rela- 
tively harder, but the source detection in the hard 
X-ray band above 125 keV is marginal (Fig. 13b). 
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6. Summary and Future Development 

We have successfully carried out the first high al- 
titude (40 km) balloon flight of the wide-field hard 
X-ray telescope, ProtoEXISTl, employing the first 
generation of an advanced CZT imager. The CZT 
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ment (256 cm 2 with 2.5 mm pixels). The X-ray tele- 
scope performed excellently throughout the 7.5 hr 
flight. Despite a few problems in the new aspect and 
pointing system, we were able to detect Cyg X-l at 
7.2a from about an hour observation, ~ 50 min of 
which was with relatively poor aspect and so with 
smearing in the elevation axis direction. 

Encouraged by this success, we have begun the 
next phase of the program, ProtoEXIST2, which 
will have 4x finer detector spatial resolution (pixel 
pitch: 0.6 mm) and also both lower energy thresh- 
old (~ 5 - 8 keV) and better energy resolution (~ 
2 keV). To accomplish this, a new ASIC, the Nu- 
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